The influence of the optical Stark effect on chiral tunnehng in graphene 
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The influences of intense coherent laser fields on the transport properties of a single layer graphene 
are investigated by solving the time-dependent Dirac equation numerically. Under an intense laser 
field, the valence band and conduction band states mix via the optical Stark effect. The chiral 
symmetry of Dirac electrons is broken and the perfect chiral tunneling is strongly suppressed. These 
properties might be useful in the fabrication of an optically controlled field-effect transistor. 
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Graphene has attracted much attention due to its re- 
markable electronic properties [l|43|- The low-energy 
quasiparticles, which have linear dispersion and nontriv- 
ial topological structure in their wave function, can be 
described by using a Dirac-like equation. This unique 
band structure of graphene leads to many important po- 
tential applications in nanoelectronics [^-|9|. 

One of the peculiar transport phenomena in graphene 
is the chiral tunneling [J, [s], |lO| . In single layer graphene 
a perfect transmission through a potential barrier in the 
normal direction is expected. This unique tunneling ef- 
fect can be explained by the chirality of the Dirac elec- 
trons within each valley, which prevents backscattering 
in general. This kind of reflectionless transmission is in- 
dependent of the strength of the potential, which lim- 
its the development of graphene-based field-effect tran- 
sistors (FET) [j|. The perfect transmission can be sup- 
pressed effectively when the chiral symmetry of the Dirac 
electrons is broken. For instance, in a magnetic field, a 
quantized transmission can be observed in graphene p- 
n Junctions [ll|. Recently, Elias et. al. proposed that 
the hydrogenation could convert the semimetal graphene 
into an insulator material |12l |. 

The intense optical field can also break the chiral sym- 
metry of Dirac electrons in graphene, e.g., Fistul and Efe- 
tov have shown that when the n-p Junctions in graphene 
is irradiated by an electromagnetic field in the resonant 
condition, the quasiparticle transmission is suppressed 
[13| . The optical field control on carrier transport of- 
fers several advantages. Optical fields can control not 
only the charge carriers but also the spin carriers, es- 
pecially which can be performed over femtosecond time 
scale. Another fundamental method of optical control is 
the optical Stark effect (OSE) [IMi. The OSE in tra- 
ditional semiconductors is due to a dynamical coupling 
of excitonic states by an intense laser field. The OSE 
have shown many useful applications in optoelectronics 
and spintronics [l9l - l23| . 

In graphene, the valence band and conduction band 
states can also mix strongly via OSE. Thus the chirality 



of Dirac electrons will be completely changed, or even 
disappear. Unlike the resonant case 13|, in OSE the 
coherent excitons are virtual excitons, which exist only 
when the optical field is present. Thus the light-induced 
shift lasts only for the duration of the pump pulse, which 
allows for optical gates that might only exist for femtosec- 
onds. Furthermore, since there is no real absorption in 
the nonresonant case, the absorption of photons is quite 
small and low power consumption is expected. 

In this Letter, we study the tunneling rate of Dirac 
electrons in graphene through a barrier with an intense 
electromagnetic field. We consider a rectangular poten- 
tial barrier with height Vq, width D in the X direction, 
and infinite length in the Y direction [see Fig. 1 (a) and 
Fig. 1 (b)]. The Fermi level (dashed lines) lies in the 
valence band in the barrier region and in the conduction 
band outside the barrier. The gray filled areas indicate 
the occupied states. The optical field is propagated per- 
pendicular to the layer surface and is linearly polarized 
along the Y direction with a detuning Aq — 2Eh — hu}. 
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FIG. 1: (color online) . (a) Schematic of the spectrum of Dirac 
electrons in single-layer graphene. The optical field is propa- 
gated perpendicular to the layer surface and and is linearly po- 
larized along the Y direction, (b) Schematic of the scattering 
of Dirac electrons by a square potential. Ba, Bi„, and Bout de- 
note the absorbing boundary, incident boundary, and output 
boundary, respectively, (c) Schematic of the one-dimensional 
Yee lattice in graphene. 



We choose Aq > to ensure that there is no interband 
absorption inside the barrier. Meanwhile, hco <C 2Ek is 
used to guarantee that the influence of the optical field 
outside the barrier can be neglected. 

Since the Coulomb interaction between electrons and 
holes in OSE is negligible when the detuning is large 



171 [19| , we did not take into account the electron- hole 
Coulomb interaction or many body effect in our calcu- 
lation. Thus, neglecting the scattering between different 
valleys, the scattering process of Dirac electrons in K 
point is described by the time-dependent Dirac equation 

^ * (r, t) = [Ho + Vo (r) I + H^*] * (r, t) , (1) 



where ^ (r, t) 
Hq = —ihvFO 



= [Cyi(r,i),CB(r, i)] is the wave function, 
' V is the unperturbed Dirac Hamiltonian, 



f = {<^x,<7y) are the Pauli matrices, vp w lOrn/s is 
the Fermi velocity, Vo(r) is the height of the potential 
barrier, I is the unit matrix, and Hint is the interaction 
Hamiltonian. H,v,t can write as 
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H, 



= -hevp [A{x, t)ax + A{y, t)ay\ = hi 
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lit) 
(2) 
where e is the electron charge and [A{x , t) , A{y , t)] = 
[^a;e*'^*,^ye*"*] are the vector potentials of the electro- 
magnetic field. When the Dirac electrons is incident on 
the barrier perdenicularly, we can rewrite Eq. ([T|) as a 
set of partial differential equations 



idCA{x,t)/dt — —ivFdCB{x,t)/dx + VoCA{x,t) 

+ Vi2{t)CB{x,t), (3) 

idCBix,t)/dt^ -ivFdCA{x,t)/dx + VoCB{x,t) 

+ V2l{t)CA{x,t). (4) 

Since the tunneling time is sub-picosecond and the po- 
tential Vi2{t) and V2i{t) vary as fast as the frequency of 
incident light beams, this scattering process is strongly 
time-dependent. In order to study such a strongly 
time-dependent scattering process, we employ the finite- 
difference time-domain (FDTD) method to solve Eq. ([3]) 
and Eq. (|4]) numerically in the time-domain 25] . In the 
traditional FDTD method, the Maxwell's equations are 
discretized by using central-difference approximations of 
the space and time partial derivatives. As a time-domain 
technique, the FDTD method can demonstrate the prop- 
agation of electromagnetic fields through a model in real 
time. Similar to the discretization of Maxwell's equa- 
tions in FDTD, we denote a grid point of the space and 
time as {i,k) = {iAx,kAt) [see Fig. 1(c)], and for the 
any function of space and time F{iAx,kAt) = F'^{i). 
the first order in time or space partial differential can be 
expressed as 
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Thus the Eq. ^ and Eq. ((4]) can be replaced by a finite 
set of finite differential equations 
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For computational stability, the space increment Aa:: and 
'the time increment At need to satisfy the relation Aa; > 
vpAt 25]. Furthermore, the space increment Aa; must 
far smaller than the wavelength of electrons Aa; < Ae/8, 
and the time increment At must be far smaller than the 
period of the electromagnetic field Ti . 

At the boundary Ba, one-dimensional Mur absorbing 
boundary conditions are used [26|. At the input bound- 
ary Bin, a Gaussian electronic wave packet is injected 



Ca = Cb = —^ cxp 



47r(t - to)2 



(8) 



where to and r denote the peak position and the pulse 
width, respectively. 

Thus, by solving Eq. ([7a|) and Eq. (l7b)) directly in 
the time domain we can demonstrate the propagation of 
a wave packet through a barrier in real time. Numeri- 
cal simulations are shown in Fig. 2. The following pa- 
rameters are used in our calculation: the peak position 
to = 1.5 ps, the pulse width t = 1.0 ps, the space incre- 
ment Aa; = 0.1 nm, the time increment At = 5 x 10~^ 
ps, and the height of the potential barrier Vq = 400 meV. 
When there is no pump beams, a perfect chiral tunnel- 
ing can be found [see Fig. 2(a)]. This result is consistent 
with that of Geim et. al. [4|. But when the sample is 
irradiated by an intense nonresonant laser beam, a re- 
fiected wave packet appears [see Fig. 2(d)]. The perfect 
transmission is suppressed. By analyzing the the trans- 
mitted wave packet and the refiected wave packet, we can 
obtain the tunneling rate. 

To explain the suppression of chiral tunneling. We first 
investigate the OSE in the barrier within a rotating- 
wave approximation 15|, [22, l23| . Figure 2(a) shows the 
renormalized band as a function of momentum k with 
intensity I^ = 30 MW/cm^. In the case of nonreso- 
nant excitation, huj < 2Ei, and the dressed states are 
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FIG. 2; (color online), (a) numerical simulations of a wave 
packet tunneling through a barrier without pump beams, (b)- 
(d) Time sequence of a wave packet tunneling through a bar- 
rier with pump intensity lu = 3 MW/cm^, Ao = 5meV, and 
D = 300 nm. The light grey shows the barrier area. 
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graphene can be broken. For instance, at very small de- 
tuning, the wave functions of these dressed states can 
be approximately written as the superposition of un- 
perturbed conduction and valence wave function, vj/ = 
('f+ + *_)/a/2 = (1,0). These dressed states are not 
the eigenstates of the helicity operator. The chiral sym- 
metry is broken and perfect chiral tunneling is strongly 
suppressed. Numerical results are shown in Fig. 2(c) 
with pump intensity /„ = 30 MW/cm^ and D = 300 
nm. From Fig. 2(c) we can find that the transmis- 
sion is strongly suppressed, even with lager detuning 
(e.g., Aq = 10 meV, the transmittance is about 0.025). 
When detuning increases, the light-induced mixing be- 
comes weak [see Fig. 2(b)], the reflectance decreases, 
and the transmittance increases. Fig. 2(d) shows the 
transmittance as a function of pump intensity with dif- 
ferent barrier widths. The strong laser field can enhance 
band mixing and reduce the transmittance. From Fig. 
2(d) we also see that the wide barrier can prolong the 
interaction time between electrons and photons, reduce 
the tunneling rate, and lower the threshold of the pump 
laser power. 

In conclusion, we have calculated the influence of the 
OSE on the chiral tunneling in graphene by using the 
FDTD method. We find that perfect tunneling can be 
strongly suppressed by the optically induced band mix- 
ing, even at large detuning. These properties might be 
useful in device applications, such as the fabrication of 
an optically controlled field-effect transistor that has ul- 
trafast switching times and low power consumption. 

This work was supported by the NSFC Grant Nos. 
10904059, 10904097, and 11004199, the NSF from Jiangxi 
Province 2009GQW0017, the Open Research Fund of 
State Key Laboratory of Milhmeter Waves No. K200901. 



FIG. 3: (color online). (a) Sketch of the renormalized 
band energies (solid lines) and the unperturbed band energies 
(dashed lines) as a function of momentum k. (b) Sketch of 
the fermion distribution function Uk as a function of momen- 
tum k. (c) The reflectance (red circles) and the transmittance 
(black squares) of the barrier as a function of the detuning for 
/„ = 30 MW/cm^ and D = 300 nm. (d) The transmittance 
as a function of pump intensity for Ao = 5 meV with different 
barrier width. 



blue shifted. With increasing detuning, the light-induced 
shift decreases, and the dressed states asymptotically ap- 
proach the unperturbed states. The intense electromag- 
netic field can also induce a strong band mixing. Near the 
absorption edge, a maximum fermion distribution func- 
tion rifc sa 0.44 can be observed [see Fig. 1(b)]. 

Under intense light beams, the dressed states are 
strongly mixed with valence states and conduction states. 
Therefore, the chiral symmetry of Dirac electrons in 
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